ToMUO®!] 3 © 001] 
PLANT SCIENCE 



ORIGINAL RESEARCH ARTICLE 

published: 11 June 2013 
doi: 10.3389/fp!s.2013. 00189 




Identification of novel transcription factors regulating 
secondary cell wall formation in Arabidopsis 

Hua Cassan-Wang , Nadia Goue, Mohammed N. Saidi f , Sylvain Legay t , Pierre Sivadon t , 
Deborah Goffner* and Jacqueline Grima-Pettenati* 

Laboratoire de Recherche en Sciences Vegetales, UMR5546, Centre National de la Recherche Scientifique, Universite Toulouse III, UPS, Toulouse, France 



Edited by: 

Maurice Bosch, Aberystwyth 
University, UK 

Reviewed by: 

Akiyoshi Kawaoka, Nippon Paper 
Industries Co., Japan 
Christian Dubos, Institut National de 
la Recherche Agronomique, France 

*Correspondence: 

Jacqueline Grima-Pettenati, 
Laboratoire de Recherche en 
Sciences Vegetales, UMR5546, 
Centre National de la Recherche 
Scientifique, Universite Toulouse III, 
UPS, 24 Chemin de Borde Rouge, 
BP 42617 Auzeville, 31326 
Castanet-Tolosan, Toulouse, France 
e-mail: grima@lrsv. ups-tlse. fr 
t Present Address: 
Mohammed N. Saidi, Laboratoire 
des Biotechnologies Vegetales 
Appliquees a /Amelioration des 
Cultures, Ecole Nationale 
d'lngenieurs de Sfax, Sfax, Tunisie; 
Sylvain Legay, IJPB, UMR1318 
INRA-AgroParisTech, Versailles, 
France; 

Pierre Sivadon, UMR 5254 IPREM, 
CNRS, Universite de Pau et des 
Pays de TAdour, Pau, France; 
Deborah Goffner, UMI 3189 Faculte 
de Medecine secteur Nord, 
Marseille, France 



The presence of lignin in secondary cell walls (SCW) is a major factor preventing hydrolytic 
enzymes from gaining access to cellulose, thereby limiting the saccharification potential of 
plant biomass. To understand how lignification is regulated is a prerequisite for selecting 
plant biomass better adapted to bioethanol production. Because transcriptional regulation 
is a major mechanism controlling the expression of genes involved in lignin biosynthesis, 
our aim was to identify novel transcription factors (TFs) dictating lignin profiles in the model 
plant Arabidopsis. To this end, we have developed a post-genomic approach by combining 
four independent in-house SCW-related transcriptome datasets obtained from (1 ) the fiber 
cell wall-deficient wat1 Arabidopsis mutant, (2) Arabidopsis lines over-expressing either 
the master regulatory activator EgMYB2 or (3) the repressor EgMYBl and finally (4) 
Arabidopsis orthologs of Eucalyptus xylem-expressed genes. This allowed us to identify 
502 up- or down-regulated TFs. We preferentially selected those present in more than 
one dataset and further analyzed their in silico expression patterns as an additional 
selection criteria. This selection process led to 80 candidates. Notably, 16 of them were 
already proven to regulate SCW formation, thereby validating the overall strategy. Then, 
we phenotyped 43 corresponding mutant lines focusing on histological observations 
of xylem and interfascicular fibers. This phenotypic screen revealed six mutant lines 
exhibiting altered lignification patterns. Two of them [Bel-like HomeoBox6 (blh6) and a 
zinc finger TF] presented hypolignified SCW. Three others (myb52, myb-like TF, hb5) 
showed hyperlignified SCW whereas the last one {hb15) showed ectopic lignification. In 
addition, our meta-analyses highlighted a reservoir of new potential regulators adding to 
the gene network regulating SCW but also opening new avenues to ultimately improve 
SCW composition for biofuel production. 
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INTRODUCTION 

Plant cells are enclosed in cell walls, which provide them with 
structural support and regulate growth and differentiation. There 
are two main types of cell walls: primary cell walls and secondary 
cell walls (SCWs). Primary cell walls are formed in all plant cells 
and are composed mainly of cellulose, hemicellulose, and pectin. 
SCWs are much thicker and are deposited in the inner side of pri- 
mary cell walls only in some highly specialized tissues and cell 
types such as xylem vessels and fiber cells. Lignified SCW are the 
most abundant source of renewable biomass on earth, and are 
widely used for construction, paper, and energy. In the context of 
the energetic crisis, lignocellulosic biomass has received growing 
attention as raw material for the production of second-generation 
biofuels. 

SCWs are composed of cellulose, hemicelluloses, and lignin. 
The impregnation with lignin renders SCWs waterproof and 
resistant, allowing water conduction through xylem vessels as well 



as mechanical support. On the other hand, for lignocellulosic 
biofuel production, lignin is a major negative factor prevent- 
ing hydrolytic enzymes from gaining access to cellulose and, as 
a result limits the saccharification potential. The biosynthetic 
pathways leading to SCWs formation are highly regulated at the 
transcriptional level. Tremendous progress has been made during 
the last decade supporting the existence of a complex hierar- 
chical regulatory network of transcription factors (TFs). Most 
of those belong to two large TF families: R2R3-MYB and NAC 
(NAM/ATAF/CUC) (Demura and Fukuda, 2007; Zhong and Ye, 
2007; Grima-Pettenati et al., 2012; Wang and Dixon, 2012; Zhong 
et al., 2012). Some members of the NAC TF family are key regula- 
tors of SCW formation in fibers and/or in vessels. This particular 
SCW related subgroup of NACs include the NAC SECONDARY 
WALL THICKENING PROMOTING FACTOR 1 (NST1), 
SECONDARY WALL- ASSOCIATED NAC DOMAIN PROTEIN 1 
(SND1/NST3), NST2, and the VASCULAR-RELATED NAC 
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DOMAIN (VND6 and VDN7) (Kubo et al, 2005; Mitsuda et al., 
2005, 2007; Zhong et al., 2006, 2012; Ko et al, 2007). Over expres- 
sion of any of these NACs led to ectopic lignification in cells 
that normally contain only primary cell walls (for review, see 
Grima-Pettenati et al., 2012). A double mutation of SND1 and 
NST1 resulted in loss of SCW in fibers, whereas the simultane- 
ous repression of VND6 and VND7 led to a defect in vessel SCW 
thickenings (Kubo et al., 2005; Zhong et al., 2007b). In the reg- 
ulatory hierarchical network SDN1, NST1/2, and VND6/7 are 
first-level master switches controlling downstream TF regulators 
(Zhong and Ye, 2007; Wang and Dixon, 2012; Zhong et al, 2012). 
The second layer of regulators includes many MYB TFs (MYB20, 
MYB42, MYB43, MYB46, MYB52, MYB54, MYB58, MYB69, 
MYB61, MYB63, MYB83, MYB85, and MYB103) (Zhong et al., 
2008; Ko et al., 2009; McCarthy et al, 2009; Zhou et al., 2009; 
Romano et al, 2012) as well as several other TFs like SND2, 
SND3, KNAT7, AtC3H14 Zinc finger TF (Zhong et al, 2008; 
Ko et al., 2009). Some of these are also master regulators since 
they control the biosynthesis of the three main components of 
SCW i.e., cellulose, xylan, and lignin. The discovery of this multi- 
leveled hierarchical regulatory network has been a breakthrough 
in our understanding of the regulation of the lignified SCW, 
although it is far from being complete. For instance, only a few 
TFs characterized hitherto are regulating specifically one of the 
SCW components, although three MYBs (MYB85, MYB58, and 
MYB63) were reported to be lignin-specific. In addition, our 
knowledge of the molecular mechanisms determining the het- 
erogeneous SCW deposition in different cell types, as well as 
those governing the various patterns of SCW deposition is still 
very poor. More efforts are needed to get a comprehensive pic- 
ture of the transcriptional regulation of the SCWs both from a 
fundamental and an applied perspective. 

As a step toward this goal, we searched for novel TFs poten- 
tially implicated in the control of lignin deposition. To do this, 
we set up a post-genomic approach combining four original 
in-house SCW-related transcriptomic data sets that enabled us 
to identify 80 candidates belonging to major plant TFs families 
(i.e., NAC, MYB, bHLH, Zinc finger, HomeoBox, and AP2/ERF). 
Most of them have not yet been functionally characterized. 
Histochemical analyses of the corresponding mutants revealed six 
strong candidates regulating the biosynthesis of lignin and/or the 
whole SCW biosynthetic program: BLH6 (Bel-like HomeoBox6; 
AT4G34610), HB5 (AT5G65310), HB15 (AT1G52150), MYB-like 
TF (AT3G11280), MYB52 (AT1G17950), and Zinc finger TF 
(AT3G46620). 

RESULTS 

A POST-GENOMIC APPROACH TO IDENTIFY NOVEL REGULATORY 
GENES INVOLVED IN SCW FORMATION 

In order to identify novel regulatory genes involved in SCW 
formation, we took advantage of four large scale in house tran- 
scriptomic data sets and developed a post-genomic approach. A 
flow chart of the main steps of this original strategy is described 
in Figure 1. The first SCW-related in house transcriptome dataset 
came from the Arabidopsis mutant, watl (walls are thin 1 ), which 
has little to no SCW in fibers (Ranocha et al., 2010). In this 
mutant, the transcript levels of many genes associated with the 



regulation and/or the biosynthesis of the different SCW wall 
polymers were dramatically reduced in keeping with the mutant 
phenotype. Within the genes up/or down-regulated in the mutant 
background, we identified 97 TFs including some well-known 
SCW-regulating TFs such as SND1, SNT1, andMYB46 (Table SI). 
The second transcriptome dataset was comprised of 240 TFs 
(Table S2) that exhibited de-regulated expression in Arabidopsis 
lines over-expressed the SCW-master activator EgMYB2, which 
is a Eucalyptus R2R3 MYB TF highly expressing in xylem cells 
undergoing SCW thickening (Goicoechea et al., 2005). EgMYB2 
is able to activate the promoters of lignin (Goicoechea et al., 
2005), cellulose, and xylan biosynthetic genes (Zhong and Ye, 
2009), leading to thicker SCW in EgMYB2 over-expressing 
lines in tobacco (Goicoechea et al., 2005; De Micco et al, 
2012). Moreover, the closest orthologs of EgMYB2 in Arabidopsis 
AtMYB46 and AtMYB83 encode for master regulators capable of 
activating the whole SCW biosynthetic program (Zhong et al., 
2007a; McCarthy et al, 2009), and EgMYB2 was able to comple- 
ment the myb46-myb83 double mutant (Zhong et al., 2010). The 
third transcriptome dataset included 309 TFs (Table S3) dereg- 
ulated in Arabidopsis lines over-expressing the SCW-repressor, 
EgMYBl (Legay et al., 2010). EgMYBl over-expressors exhibited 
fewer lignified fibers particularly in the interfascicular zones and 
reduced SCW thickenings. Klason lignin content was moderately 
but significantly reduced and decreased transcript accumulation 
was observed for genes involved in the biosynthesis of lignins, cel- 
lulose, and xylan (Legay et al., 2010). Finally, the fourth dataset 
was composed of 87 TFs (Table S4) that were the Arabidopsis 
orthologs of Eucalyptus TFs preferentially expressed in differen- 
tiating xylem (Rengel et al., 2009), a tissue that is particularly rich 
in cells undergoing SCW deposition and lignification. Altogether, 
these four transcriptomic datasets allowed us to identify a total 
of 502 candidate TFs. To narrow down the number of candidates 
for functional validation, we selected 186 that were identified in 
two datasets (Table S5). It should be noted that 43 of those were 
found in three data sets and only three were common to the four 
datasets bHLH5 (AT5G46760), IAA9 (AT5G65670), and AP2 TF 
RAP2.2 (AT3G14230). 

CROSS-COMPARISON WITH PUBLICLY AVAILABLE MICR0ARRAY DATA 

To further narrow down the selection of the 186 candidate TFs for 
further functional analysis, we examined their in silico expression 
patterns using Genevestigator (Hruz et al., 2008). We restricted 
our list to genes that were preferentially and/or highly expressed 
in situations in which SCW formation is prevalent i.e., in xylem, 
the basal part of the inflorescence stem, and/or in cell suspension 
cultures undergoing in vitro SCW formation (Kubo et al., 2005). 
This in silico expression screen allowed us to obtain a final list of 
80 candidate SCW TFs (Table 1). 

The most represented families in decreasing order were the 
MYB TF family (containing 19 members), the AP2 ERF (Ethylene 
Response Factor) TF family (12), the HB (HomeoBox protein) TF 
family (9), and the Aux/IAA family (7) (Figure 2). The other TFs 
belonged to the WRKY (5), NAC (5), bHLH (5), Zinc finger TF 
(5), bZip (3), and ARF (2) families, respectively. 

It is noteworthy that 16 of the 80 candidates were already 
shown to regulate SCW formation. They include eight MYB 
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Four in-house transcriptomic datasets: 502 TFs 




1. Arabidopsis mutant wat1 

2. Arabidopsis line overexpressing EgMyb2 

3. Arabidopsis line overexpressing EgMyM 

4. Aradidopsis orthologs of Eucalyptus xylem expressed TFs 



186 TFs 

Expressed in more than one dataset 




Expressed in xylem 



in vitro TE induction 



Genevestigator 



80 selected candidate transcription factors 
(of which 16 known SCW regulators) 



43 T-DNA / RNAi lines phenotyped 



6 mutants with altered SCW 



FIGURE 1 | Overall strategy to identify transcription factors (TFs) 
involved in SCW formation. Four in-house SCW formation-related 
transcriptomic datasets were crossed to select TFs present in more 
than one dataset/experimental condition. These TFs were further 
screened against publicly available large-scale transcriptomic datasets, 
to select those highly or preferentially expressed in the organs and/or 
tissues of interest. This led to a list of 80 candidate genes, of which 



we phenotyped the 42 T-DNA insertion mutants and/or RNAi 
transgenic plants available. The four in house SCW related 
transcriptomic datasets include (1) the fiber SCW-deficient watl 
Arabidopsis mutant, (2) Arabidopsis lines over-expressing the SCW 
master activator EgMYB2, (3) Arabidopsis lines over-expressing the 
SCW master repressor EgMYBI, and (4) Arabidopsis orthologs of 
Eucalyptus xylogenesis-related genes. 



TF (MYB20, MYB42, MYB46, MYB52, MYB61, MYB63, MYB75, 
MYB85, and MYB103) (Zhong et al, 2007a, 2008; Ko et al, 2009; 
Zhou et al., 2009; Ohman et al, 2013), four NAC TF (SND1, 
NST1 XND1, and VNI2) (Mitsuda et al., 2005; Zhong et al., 2006; 
Zhao et al., 2008; Yamaguchi et al., 2010), three Homeodomain 
containing TF (HB14, HB15, and KNAT7) (McConnell et al., 
2001; Emery et al., 2003; Kim et al., 2005; Zhong et al, 2008) 
and one WRKY TF (WRKY12) (Wang et al, 2010) (complete 
list in Table 1). This significant proportion of characterized SCW 
related TF in our final list validates well the strategy used in this 
study. For example, MYB46 (Zhong et al, 2007a), MYB85 (Zhong 
et al, 2008), and WRKY12 (Wang et al., 2010) were present in 
three of the four transcriptomic datasets. 

PHEN0TYPES OF TF T-DNA MUTANT OR RNAi TRANSGENIC LINES 

We then collected and characterized publicly available T-DNA 
mutant lines or RNAi transgenic lines that corresponded to 43 of 
the 80 candidate genes (Table S6). The information concerning 



the different lines including T-DNA insertion position and in- 
house databases source is presented in Table S6. Phenotyping 
was performed on 20 cm-high mutant stems grown in short- 
day growth conditions. Under these conditions, the basal part 
of the stem abundantly develops cells undergoing SCW thicken- 
ing (xylem vessel cells, xylary fiber cells, and interfascicular fiber 
cells). Histological analyses of SCW were performed using the 
natural auto fluorescence of phenolic compounds under UV-light 
as well as phloroglucinol-HCl staining, which is indicative of the 
lignin content. We found significant alteration of lignin profiles 
in six mutant lines corresponding to two MYB TFs: MYB like 
TF (AT3G11280) and MYB52 (AT1G17950), three HomeoBox TF 
HB5 (AT5G65310), BLH6 (AT4G34610), and HB15 (AT1G52150) 
and a Zinc finger TF (AT3G46620), although the overall organiza- 
tion of vascular bundles and interfascicular fibers was not altered 
in these six mutant lines (Figures 3, 4). 

Under UV-light the intensity of auto-fluorescence was lower in 
zinc finger TF (Figure 3B) and in blh6 (Figure 3C) mutant lines in 



www.frontiersin.org 



June 2013 | Volume 4 | Article 189 | 3 



Cassan-Wang et al. 



Novel regulators of lignified secondary walls 



Table 1 | List of 80 candidate TFs obtained after cross-comparison of transcriptomic datasets. 



ID Annotation in house data sets (Log2Ratio) Publicly available data No. of ordered 
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Table 1 | Continued 
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in house data sets (Log Ratio) 


Publicly available data 


No. of ordered 
mutants 


wat1 


ut 

tgivf idz 


Ut 

tgivi Y ts I 


tg 

Xylem 


TE induction 6d* 


Xylem** 


A I 3b4o13(J 


MYd4o 


3.0 


1.4 


1.8 




0.1 


High_s 


0 


ATI Cj 17950 


iuivdo 

MYB5Z 


-2.1 


2.3 


2.4 




3.3 


Medium_s 


2 


atcp cm on 
Al buoy /oU 


IVI Y bby 


4.7 


1.0 


0.8 




2.7 


High_p 


1 


Al luUyo4U 


IVI Y bo 1 


-2.4 




1.8 




no probe 


no probe 


0 


Al lu /»loU 


IVI Y boo 


-5.9 


1.3 






2.3 


Medium_s 


0 


ATz.36650 


RJIVD7E 

MYB75 




-1.3 


-1.9 




0.1 


Low_u 


0 


Al 3bbUUoU 


MYd// 


-3.2 


1.4 


0.8 




-0.5 


Medium_u 


0 


AT4Cj22680 


MYB85 


-2.1 


2.2 


1.6 




2.1 


Low_u 


1 


Al bbUb/yU 


myb-like TF 


-1.9 


0.7 


0.9 




1.7 


High_s 


0 


a tit r 1 1 Tonn 

Al bbl /30U 


myb-like TF 


-1.8 


0.8 






1.2 


Medium_u 


0 


Al 3u11zoU 


myb-like TF 


1.9 




0.6 




4.5 


High_p 


3 


ati 1 nrn 

Al 1 ulzzbU 


NALUU/, VNU4 


-1.9 








5.0 


Low_u 


1 


ati r'O'mn 
Al I u Jz/ /U 


MAPA11 Cl\ir\1 MCTO 

IMACU Iz, oIMLm, 1Mb 1 o 


-2.6 




1.8 




0.0 


Medium_s 


0 


Al Zu4o/ /U 


l\l APA/IO |\|CT1 

IMAUU4J, No 1 1 


-0.8 






X 


no probe 


no probe 


0 


ATbCa 64530 


NAC104, XIMD1 






0.9 


X 


3.6 


High_s 


0 


ate pioi on 
Al buloloU 


IMAUoo, VNU-IN 1 tKAU 1 IIMuZ ( V IMIz) 




-0.7 




X 


3.3 


High_u 


0 


Al 4bo / /bU 


AIMTCPI IP\/1CMTA /AMU 

AIN 1 LbUIVILIM IA (AN 1 ) 


-5.4 




1.0 




0.1 


Low_p 


0 


AT1 POl /icn 

Al 1 bz l4b(J 


scarecrow-iiKe ir 1 (oLLi) 




1.2 




X 


2.2 


High_p 


1 


Al Zb4/U/U 


squamosa TF-like 1 (SPL1) 




0.8 




X 


0.7 


High_u 


0 


Al zu44/4b 


1A/DI/V10 


-2.9 


2.3 


2.7 




-0.1 


High_s 


1 


Al zuoUZbU 


VvnKYzb 


2.6 


-0.7 


-0.6 




-1.2 


Medium_u 


2 


Al 4bzobbU 


VvnKYzy 


-7.7 


2.1 


1.6 




no probe 


no probe 


1 


AT1 HQHQA fl 

A I luoUo4U 


Vvnl\Y4U 


-3.0 






X 


-0.4 


High_u 


I 


AT3G01970 


WRKY45 


3.4 


-0.8 


-1.2 




-0.5 


Medium_u 


0 


AT5G 13080 


WRKY75 


3.5 


-0.8 


-0.9 




-0.3 


Low_ns 


0 


AT2G28510 


zinc finger Dof-type 


-1.7 


-0.8 


-0.9 




0.2 


Low_ns 


0 


AT1G68360 


zinc finger TF 


-3.7 


1.8 


2.0 




1.1 


High_s 


0 


AT2G40140 


zinc finger TF (CCCH-type) 


-1.8 


0.6 






0.3 


High_u 


0 


AT3G 55980 


zinc finger TF (CCCH-type) 


-3.6 






X 


0.5 


High_u 


2 


AT3G46620 


zinc fingerTF-n129 


-3.4 








1.6 


High_u 


2 



*TE induction 6d, 6 days after induction (1 uM brassinolide and 10 mM H3B03) when tracheary element were actively forming. 

**Xylem in hypocotyl of adult plant; High, Medium, and Low for their expression level; _s, specific expression in xylem; _p, preferential expression; _u, ubiquitous 
expression; _ns, higher expression in non-xylem cells. Bold indicates known genotypes, and corresponding secondary wall or xylem cell identity phenotypes. 



both vascular bundles and interfascicular regions as compared to 
the control (Figure 3A), suggesting a global decrease in phenolic 
compound deposition. The SCW in xylem vessels, xylary fibers, 
and interfascicular fibers were observed in more detail using 
phloroglucinol-HCl staining. Little to no SCW was deposited in 
xylary fibers (Figures 3H,I) and moreover, SCW thickness was 
largely reduced in interfascicular fibers (thin and weak phloroglu- 
cinol staining) as compared to wild-type (Figure 3G) suggesting 
that these lines were hypolignified. 

Auto-fluorescence under UV light was more intense in myb 
like TF (Figure 3D), hb5 (Figure 3E), and myb52 (Figure 3F) 
lines than in controls (Figure 3A), especially in the interfas- 
cicular region, suggesting an increased deposition of phenolic 
compounds and possibly lignins. This was further confirmed 
by a massive and intense phloroglucinol staining indicating an 
enhanced lignin deposition in the interfascicular fiber and xylary 



fiber cells of these mutants (Figures 3J-L). Extra-layers of cells 
with lignified SCW were detected in the external layers of both 
interfascicular fibers and metaxylem vessels in two lines myb- 
like TF (Figure 3J, green arrows) and myb52 (Figure 3L, green 
arrows) as compared to the control (Figure 3A). This observa- 
tion suggests that secondary xylem formation was enhanced and 
appeared earlier than in wild-type. A strong fluorescent signal 
was also detected in the phloem cap cells (Figures 3D-F, blue 
arrows) in all three highly auto-fluorescent lines suggesting a 
transition of phloem cap cells to phloem sclereids (highly lig- 
nified) which was further confirmed by phloroglucinol staining 
(Figures 3J-L, blue arrows). Similarly, auto -fluorescent signals 
(Figure 3F) and strong phloroglucinol-HCl staining (Figure 3L, 
pink arrow) were detected in the epidermal cells of some 
myb52 T-DNA insertion lines revealing an ectopic deposition of 
lignin. 
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15 



20 




FIGURE 2 | Distribution of the 80 SCW candidate TFs according to 
family. Amongst them, four families have been shown previously to 
regulate SCW formation: MYB, NAC, HB (Homeodomain containing TF), 
and WRKY; three TF families are involved in hormone signaling including 
ARF Aux/IAA, and AP2 ERF TF families. "Others" include the following TF 
families: LIM, CCAAT binding, HMG (High Mobility Group, belonging to a 
transcription complex), AINTEGUMENTA (ANT; AP2 like TF), Scarecrow-like 
TF1 (SCL1), and Squamosa-likel (SPL1). 



Both auto-fluorescence and phloroglucinol staining of stem 
sections of hbl5 (Figures 4B,D) showed that large parenchyma 
cells adjacent to the inner side of the interfascicular fibers (red 
arrow), as well as smaller xylem parenchyma cells surrounding 
the protoxylem (yellow arrow) exhibited lignified SCW. The cor- 
responding cells in the control have non-lignified primary walls 
(Figures 4A,C). As compared to the control, extra layers of cells 
with lignified SCW were present in the most external rows of 
the interfascicular fibers and xylem (Figures 4B,D, green arrows) 
suggesting an enhanced and early formation of secondary xylem. 
Moreover, both xylary and interfascicular fibers in hbl5 lines 
exhibited both a more intense auto-fluorescence and staining by 
phloroglucinol than that of the control suggesting higher lignin 
content. 

The overall growth behavior of the mutants did not differ 
significantly from the controls, except the bolting and flowering 
time were altered in three of the mutant lines. Hypolignified blh6 
and zinc finger lines bolted and flowered earlier than controls 
(Figures 5A,C) whereas the hyperlignified hbl5 line exhibited 
delayed bolting and flowering (Figures 5B,C). In addition, hbl5 
mutants exhibited aerial rosettes at the base of the lateral inflores- 
cence branches instead of growing cauline leaves as in wild-type 
plants (Figure 5D). 

CO-EXPRESSION ANALYSIS OF CANDIDATE TFs GENES 

Since it is known that transcriptionally coordinated genes tend 
to be functionally related (Ruprecht and Persson, 2012), we per- 
formed co-expression analyses for the six candidate genes in order 
to further validate their role in controlling SCW synthesis and 
get some clues about their function. The co-expression genes lists 
were generated using the Genevestigator platform (https://www. 
genevestigator.com), Arabidopsis co-expression data mining tools 
and GeneCAT. All six candidate TFs were co-expressed with genes 
related to cell wall formation (Tables 2, 3 and Tables S7-S10) 



albeit to different extents ranging from 10 to 66% of SCW asso- 
ciated genes amongst the 50 first co-expressed genes. The most 
remarkably high co-expression profiles were found for MYB52 
and BLH6. 

Among the top 50 genes co-expressed with BLH6, 30 (60%) 
were related to SCW formation (Table 2). Notably, they included 
genes involved in the biosynthesis of the three main cell 
wall polymers i.e., three main SCW cellulose synthases genes 
(CESAs) (IRX5/CESA4, IRX3/CESA7, and IRX1/CESA8), three 
genes implicated in xylan biosynthesis [IRX10/GUT2, GUX1 
(Glucuronic acid substitution of xylanl), and XTH32 (xyloglu- 
can endotransglucosylase/hydrolase 32) (Brown et al, 2009; 
MacMillan et al, 2010)] as well as two laccase genes {IRX12/LAC4 
and LAC17) involved in monolignol polymerization (Berthet 
et al., 2011). Moreover, several SCW master transcriptional regu- 
lators were also co-expressed with BLH6, such as MYB46, MYB83, 
and SND2 (Zhong et al., 2008; Zhong and Ye, 2012). 

The co-expression analysis of MYB52 revealed that 33 of the 
top 50 genes (66%) co-expressed with MYB52 were related to 
SCW formation. The list of co-expressed genes included genes 
encoding key biosynthetic enzymes of the main SCW poly- 
mers i.e., cellulose (IRX3/CESA7, IRX1/CES8, IRX6/COBL4, and 
CTL2), hemicelluloses (IRX10/GUT1, IRX8, IRX9, IRX15, and 
IRX15-L), and possibly lignins (LAC2 and LAC10). MYB52 was 
also co-expressed with many SCW regulators such as MYB54, 
MYB85, MYB69, MYB 103, SND2/NAC073, and XND1/NAC104 
(Xylem NAC Domainl). 

DISCUSSION 

Functional genomics approaches developed during the last 
decade have generated numerous candidate genes related to SCW 
formation in Arabidopsis and other plant species. Whereas these 
large individual gene lists make difficult the choice of the most 
promising candidates for the further functional validation, meta- 
analyses combining multiple transcriptomic data sets offer a new 
way to reveal some core regulators. 

By cross-comparing four SCW-related transcriptomic 
datasets, we selected 186 TFs present in at least two experimen- 
tal conditions. Since these datasets came from very different 
backgrounds (mutants and over-expressors of SCW regulatory 
genes as well as orthologs of Eucalyptus xylem expressed genes), 
the selection of genes appearing in more than one dataset likely 
helped us to identify "core regulators" of SCW formation but 
might also have filtered out some more specific regulators. We 
further restricted the candidate gene list by including in silico 
analyses of their expression making the hypothesis that TFs 
expressed highly or preferentially in xylem tissues and/or during 
tracheary elements formation would be the most promising 
candidates. Indeed this strategy was successful since among 
the 80 genes that came out, 16 have already been reported 
to be regulators of the SCW They included, for instance, the 
master regulators SND1 and MYB46 as well as the lignin-specific 
MYB85. 

Forty-three mutant lines were phenotyped but only six exhib- 
ited a notable cell wall phenotype. This high proportion of 
mutants without phenotype is not surprising since many mutants 
targeting only one TF are known to yield mild to no phenotype 
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FIGURE 3 | Stem cross sections of five T-DNA mutant lines presenting 
either hypo or hyperlignified SCWs. Sections of wild-type plant and T-DNA 
mutants were observed under UV light (A-F) or stained with 
phloroglucinol-HCI (G-L). Phloem cap cells and ectopic lignification in 
epidermal cells are indicated by blue and pink arrows, respectively. 



Observations were made at the basal part of inflorescence stems at the 
stage of newly formed green siliques, about two weeks after bolting, when 
the inflorescence stems reached 20 cm height, if, interfascicular fiber; xf, 
xylary fiber; mx, metaxylem; px, protoxylem; sx, secondary xylem; ep, 
epidermis. Scale bar: 20|liti. 
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FIGURE 4 | Stem sections of the hb15 T-DNA mutant. Auto-fluorescence 
observed under UV light of wild type (A) and hb15 mutant (B) 
Phloroglucinol-HCI staining of lignin in wild type (C) and hb15 mutant (D). 
Ectopic lignification in large parenchyma cells and in small parenchyma cells 
surrounding protoxylem is indicated by red arrows and yellow arrows, 
respectively; precocious secondary walled secondary xylem formation is 
indicated by green arrow, if, interfascicular fiber; xf, xylary fiber; mx, 
metaxylem; px, protoxylem; sx, secondary xylem; ep, epidermis. 
Scale bar: 20 M-m. 



(Okushima et al., 2005; Overvoorde et al., 2005; Jensen et al., 
2011; Ruprecht et al., 2011). This is particularly true for multi- 
gene families TFs where functional redundancy prevents the 
observation of distinct phenotypes in knock-out mutants. This 
is indeed the case for a large proportion of the SCW regulators 
characterized hitherto including some of the sixteen highlighted 
here. For instance, whereas a single mutant of the SCW mas- 
ter transcriptional activator MYB46 did not exhibit any cell wall 
phenotype, the double knock out mutant myb46/myb83 with its 
closest ortholog MYB83 showed a severe reduction of SCW thick- 
ness (Zhong et al., 2007a; McCarthy et al., 2009). Therefore, genes 
for which the corresponding single mutants exhibited no phe- 
notype in this study may still be interesting candidates taking 
part in the regulation of SCW formation. Further experiments 
using over-expressors and/or mutants of two or more paralog 
genes would increase the probability of obtaining informative 
phenotypes and insight into their functions. Our in silico analyses 
pointed out some very promising genes which should be further 
characterized using such approaches. 




blh6 zinc hb5 
finger TF 



FIGURE 5 | Comparison of flowering time between wild-type and 
mutants plants. (A) Early flowering hypo-lignified line blh6. (B) Delayed 
flowering of hyperlignified line bh5. (C) Flowering time comparison 
between wild type and hypo or hyperlignified mutants. Significant and very 
significant statistical differences to wild type are represented by * or **, 
respectively (Student f-test). (D) Aerial rosette formation on the 
inflorescence of hb5 mutants. Scale bar: 5cm. DAG, days after 
germination. Arrows point to aerial rosettes. 



The most abundantly represented TF family in our list was the 
MYB family (19 members) of which eight (belonging to the R2R3 
subfamily) have already been shown to regulate either the phenyl- 
propanoid pathway and/or the SCW formation. It is the case for 
MYB46 (Zhong et al, 2007a), MYB63 (Zhou et al, 2009), MYB85 
(Zhong et al., 2008), and MYB103 (Ohman et al, 2013). We 
phenotyped myb52 insertion lines that exhibited a strong hyper- 
lignification phenotype, thus suggesting that MYB52 could be a 
repressor of the lignin biosynthesis and possibly of the whole 
SCW formation. This result is in apparent contradiction with a 
previous study showing that the dominant repression of MYB52 
caused a severe reduction in SCW thickening in both interfas- 
cicular fibers and xylary fibers of the inflorescence stem (Zhong 
et al, 2008). The authors concluded that MYB52 was an acti- 
vator of the SCW although no phenotype was detectable when 
over-expressed. A likely explanation to these apparent discrepan- 
cies is that MYB52 encodes a transcriptional repressor as clearly 
suggested by our knock-out mutant phenotype and therefore its 
dominant repression would result in a stronger transcriptional 
repression. MYB52 appeared to be tightly co-expressed with 
MYB54 and WAT1. It is also co-expressed with several cellulose 
and xylan biosynthetic genes and with MYB85, a specific regu- 
lator of the lignin biosynthesis (Zhong et al., 2008). Altogether, 
these results suggest for MYB52 a repressor role of the whole 
SCW program although this needs to be supported by further 
experiments. 

Besides these canonic R2R3 MYBs, four MYB-like proteins 
were present in the candidate list and one mutant was analyzed. 



Frontiers in Plant Science | Plant Biotechnology 



June 2013 | Volume 4 | Article 189 | 8 



Cassan-Wang et al 



Novel regulators of lignified secondary walls 



Table 2 | Top 50 co-expressed genes with BLH6 using Genevestigator platform. 



Probeset 


AG I 


Gene symbol 


Score 


Description 


253247_at 


AT4G34610 


BLH6 


1 


BELI-like homeodomain 6 


252025_at 


AT3G52900 




0.75 


Family of unknown function (DUF662)* 


250600_at 


AT5G07800 




0.73 


Flavin-binding monooxygenase family protein 


264573_at 


AT1G05310 




0.71 


Pectin lyase-like superfamily protein 


245657_at 


AT1G56720 




0.71 


Protein kinase superfamily protein 


251069_at 


AT5G01930 


MAN6 


0.70 


Glycosyl hydrolase superfamily protein 


263629_at 


AT2G04850 




0.70 


Auxin-responsive family protein 


250933_at 


AT5G03170 


FLA 11 


0.69 


FASCICLIN-like arabinogalactan-protein 11 


262796_at 


AT1G20850 


XCP2 


0.68 


xylem cysteine peptidase 2 


245864_at 


AT1G 58070 




0.68 


Unknown protein 


258938_at 


AT3G 10080 




0.68 


RmlC-like cupins superfamily protein 


251131_at 


AT5G01190 


LAC 10 


0.68 


laccase 10* 


246827_at 


AT5G26330 




0.67 


Cupredoxin superfamily protein* 


246439_at 


AT5G 17600 




0.67 


RING/U-box superfamily protein 


251249_at 


AT3G62160 




0.66 


HXXXD-type acyl-transferase family protein 


257757_at 


AT3G 18660 


PGSIP1.GUX1 


0.66 


Glucuronic acid substitution of xylan 1 


264184_at 


AT1G54790 




0.65 


GDSL-like Lipase/Acylhydrolase family protein* 


253876_at 


AT4G27430 


CIP7 


0.65 


COP1-interacting protein 7 


259688_at 


AT1G63120 


RBL2 


0.65 


RHOMBOID-like 2 


267094_at 


AT2G38080 


LAC4,IRX12,LMC04 


0.65 


Laccase/Diphenol oxidase family protein 


250664_at 


AT5G07080 




0.64 


HXXXD-type acyl-transferase family protein 


250322_at 


AT5G 12870 


MYB46 


0.64 


myb domain protein 46 


250120_at 


AT5G 16490 


RIC4 


0.64 


ROP-interactive CRIB motif-containing protein 4 


251297_at 


AT3G62020 


GLP10 


0.64 


germin-like protein 10 


256367_at 


AT1G66810 




0.64 


Zinc finger C-x8-C-x5-C-x3-H type family protein 


265277_at 


AT2G28410 




0.63 


Unknown protein 


246425_at 


AT5G 17420 


IRX3,MUR10,CESA7 


0.63 


Cellulose synthase family protein 


266424_at 


AT2G41330 




0.62 


Glutaredoxin family protein 


253798_at 


AT4G28500 


ANAC073,SND2 


0.62 


NAC domain containing protein 73 


260867_at 


AT1G43790 


TED6 


0.62 


Tracheary element differentiation-related 6 


251050_at 


AT5G02440 




0.61 


Unknown protein 


264493_at 


AT1G27440 


GUT1,GUT2,IRX10 


0.61 


Exostosin family protein 


261653_at 


AT1G01900 


SBTI1.1 


0.61 


Subtilase family protein 


262922_at 


AT1 G79420 




0.61 


Protein of unknown function (DUF620) 


266244_at 


AT2G27740 




0.61 


Family of unknown function (DUF662)* 


257896_at 


AT3G 16920 


CTL2 


0.60 


chitinase-like protein 2 


249070_at 


AT5G44030 


IRX5,NWS2,CESA4 


0.60 


cellulose synthase A4 


247590_at 


AT5G60720 




0.60 


Protein of unknown function, DUF547* 


267037_at 


AT2G38320 


TBL34 


0.60 


TRICHOME BIREFRINGENCE-LIKE 34* 


247648_at 


AT5G60020 


LAC 17 


0.60 


laccase 17 


257151_at 


AT3G27200 




0.60 


Cupredoxin superfamily protein* 


256155_at 


AT3G08500 


MYB83 


0.60 


myb domain protein 83 


254618_at 


AT4G 18780 


IRX1,CESA8,LEW2 


0.60 


Cellulose synthase family protein 


251630_at 


AT3G57420 




0.60 


Protein of unknown function (DUF288)** 


260914_at 


AT1G02640 


BXL2 


0.60 


beta-xylosidase 2** 


ZOOJO/ d I 


Ml OU 4JOU 


SRF7 


0 59 


^TRI IRRFI \d rprpntnr familv 7 


263841 _at 


AT2G36870 


XTH32 


0.59 


xyloglucan endotransglucosylase/hydrolase 32 


259657_at 


ATI G 55 180 


PLDEPSILON.PLDALPHA4 


0.59 


phospholipase D alpha 4 


255150_at 


AT4G08160 




0.59 


xylanase, glycosyl hydrolase family 10 protein* 


246512_at 


AT5G 15630 


IRX6.COBL4 


0.58 


COBRA-like extracellular glycosyl-phosphatidyl 
inositol-anchored protein family 


253191_at 


AT4G35350 


XCP1 


0.58 


xylem cysteine peptidase 1 


'Annotated 


involved in "xylan biosynthetic 


process" and/or "cell wall 


biogenesis " and/or 


"cell wall macromolecule metabolic process" in lair 



(http://www.arabidopsis.org). 

* 'Located in cell wall according to Tair. Bold indicates related to cell wall formation. 
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Table 3 | Top 50 co-expressed genes with MYB52 using Genevestigator platform. 



Probeset 


AGI 


Gene symbol 


Score 


Description 


255903. 


.at 


Alibi /ybu 


MYdoz,dWoz 


1 


myb domain protein 52 


245735 


at 


AT1Cj734iO 


MYB54 


0.90 


myb domain protein 54 


252550. 


.at 


ato r* a com 
Al ob4bo/U 




0.82 


nodulin MtN21 /EamA-like transporter family protein 


264559. 


at 


ATI ti 096 10 


GXM3 


0.82 


PI 1 IPI IDAM AWI A IM l\ /I CTUVI TD A MCCCD ACC O / 1 1 rcifll 

CjLUCURONOXYLAN McTHYLTRANSrcRASE 3 (DUrb79) 


248761. 


at 


Al bb4/bJb 




0.81 


Pollen Ole e 1 allergen and extensin family protein 


252211. 


at 


AT3(j50Z20 


IRX15 


0.81 


Protein of unknown function (DUF579) 


266783. 


at 


ATonoi OA 
Al ZuZal oU 


1 APO 


0.80 


laccase z 


257233. 


at 


atop 1 cncn 
Al 3ul bUbU 


IUU1U 


0.79 


lu-domain iu 


256054. 


.at 


ati f p.—; -inn 
Al 1 bl) / 1 zl) 




0.79 


Tetratricopeptide repeat (TPR)-like superfamily protein 


254277. 


at 


Ml 4uZ&uoU 


MVRQC 
IVI Y Dob 


U./o 


myb domain protein 85 


264493. 


at 


AT 1 n o~ja A n 


f2l IT1 f2l IT9 IRVin 
UUI I,UU IZ,lnA IU 


n "7"7 
u./ / 


Exostosin family protein 


250770. 


at 


ai obuwyu 


I API 1 ) 


n n 
U.// 


laccase 12* 


246425. 


.at 


ATcp 1 "7 a on 
A I bul /4zU 


idvo iv/ii iDin pccai 
lnA3,IVIUrUU I LrboA/ 


0.76 


Cellulose synthase family protein 


251009. 


at 


Al bb(Jzb4U 




0.76 


Unknown protein 


259584. 


.at 


Al I bzoUoU 




0.76 


RING finger protein 


248121. 


at 


ATCPC/ICQn 

Al b(jb4byU 


/"*AI ITIO IDVO I PTC 

uAU 1 1 z,ln ao,Lu 1 b 


0.75 


galacturonosyltransferase 12 


261928. 


at 


ati nil a on 
Al 1 bZz4oU 




n "7/i 
U. /4 


Cupredoxin superfamily protein 


254170. 


at 


AT/1 CO A A Qf"l 
A 1 4bZ44oU 




U. /4 


Rhamnogalacturonate lyase family protein 


246512. 


at 


A I bulbboU 


IKAb,UUbL4 


0.74 


UUbKA-llKe4 


264495. 


.at 


ati nioon 
Al I bz / ooU 


nlbz 


0.74 


ROP-interactive CRIB motif-containing protein 2 


245105. 


at 


Al Zb4 lb IU 




0.74 


Unknown protein 


253327. 


at 


Al 4uo34bU 


IVI Y boy 


U./4 


myb domain protein 69 


257896. 


at 


i\Top 1 coon 
A I 3U I bijzU 


PTI O 

(# 1 Lz 


0.73 


chitinase-like protein 2 


255150. 


at 


A 1 4uUo1 bU 




0.73 


glycosyl hydrolase family 10 protein* 


266714. 


at 






0.73 


No_match 


260326. 


at 


ATI (j 639 10 


IVI YB103 


0.73 


myb domain protein 103 


248907. 


at 


Ml Du40j4U 


nvVM I 


V. /z 


ppni IPPn \A/AI 1 APFTVI ATIOM 1 O aratnltrancforaca familw 
nCUULCU VVMLL MUC I Y LAI ILHM I , KJ-aCGXy ITranSTeraSe Tamily 


266488. 


at 


Al zuj4/b/U 




0.72 


Plant invertase/pectin methylesterase inhibitor family protein 


253379. 


at 


Al 4uj 6oo6\J 


DPCIDO PI IVO 

rujblr J,UU Az 


0.72 


plant glycogenin-like starch initiation protein 3 


261399. 


at 


at 1 nocin 
Al 1 uj /ybzU 




0.72 


Leucine-rich repeat protein kinase family protein* 


251478. 


at 


ATOPCQCOn 

Al jbbsbsU 


IUU13 


0.71 


lu-domain 13 


265463. 


at 


AT2G37090 


IRX9 


0.71 


GT43 glycosyl transferase43 


266708. 


at 


AT2G 03200 




0.71 


Eukaryotic aspartyl protease family protein* 


248887. 


.at 


AT5G46115 




0.71 


Unknown protein 


267414. 


at 


AT2G34790 


MEE23,EDA28 


0.70 


FAD-binding Berberine family protein 


246344. 


.at 


AT3G56730 




0.70 


Putative endonuclease or glycosyl hydrolase 


263470. 


at 


AT2G31900 


ATMY05,XIF 


0.70 


myosin-like protein XIF* 


264305. 


.at 


AT1 G78815 


LSH7 


0.70 


LIGHT SENSITIVE HYPOCOTYLS 7, Protein of unknown function 

/r~ii ice/ n\ 
(UUrb4U) 


251131. 


at 


ATcpnuon 

ai buu i iyu 


I Ann 


0.69 


laccase 10* 


251093. 


at 


ATCPni ocn 
Al buUlobU 


TDI O 

I bLJ 


0.69 


Plant protein of unknown function (DUF828)* 


247264. 


at 


atcpc/i con 
Al bUb4boU 


ANU1,AIMAU1U4 


0.69 


xylem NAC domain 1 


260430. 


.at 


Al 1 bbozUU 




0.69 


Zinc finger C-x8-C-x5-C-x3-H type family protein 


261809. 


at 


Al luUoo4U 




u.by 


Rho GTPase activating protein* 


253798. 


at 


Al 4(jzobUU 


o IM U Z , IM AO U / 3 


0.69 


NAC domain containing protein 73 


246342. 


.at 


AT3G56700 


FAR6 


0.69 


Fatty acid reductase 6 


247030. 


at 


AT5G67210 


IRX15-L 


0.69 


Protein of unknown function (DUF579) 


254618. 


at 


AT4G18780 


IRX1,CESA8,LEW2 


0.69 


Cellulose synthase family protein 


247590. 


at 


AT5G60720 




0.68 


Protein of unknown function, DUF547* 


249439. 


at 


AT5G40020 




0.68 


Pathogenesis-related thaumatin superfamily protein* 


253877. 


at 


AT4G27435 




0.68 


Protein of unknown function (DUF1218)* 


253710. 


at 


AT4G29230 


NAC075 


0.68 


NAC domain containing protein 75 



* Annotated involved in "xylan biosynthetic process" and/or "cell wall biogenesis" and/or "cell wall macromolecule metabolic process" in Tair 
(http://www.arabidopsis.org). Bold indicates related to cell wall formation. 
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Although none of these MYB-like factors has been yet reported 
as regulators of the SCW, the myb-like TF T-DNA mutant had 
a clear hyperlignification phenotype suggesting a repressor role 
of the lignin biosynthesis and/or SCW. The myb-like TF gene 
was annotated in TAIR (http://www.arabidopsis.org) as a puta- 
tive MYB domain containing TF able to interact with the gene 
product of vacuolar ATPase subunit Bl (VHA-B1). Interestingly, 
it is highly co-expressed with a newly reported gene XIP1 (XYLEM 
INTERMIXED WITH PHLOEM 1), a leucine-rich repeat receptor- 
like kinase (Table S8). The XIP1 knock-down mutants shows 
the accumulation of cells with ectopic lignification in regions of 
phloem in the vascular bundles of inflorescence stems (Bryan 
et al, 2012). 

The homeodomain containing TFs were well represented 
in the list of candidate genes with nine members. Members 
of this family have been shown to regulate procambium cell 
activities by promoting secondary walled xylem cell differen- 
tiation during vascular development. Some HD-ZIP III TF 
(HB8, PHV/HB9, PHB/HB14, REV/IFL1, and CAN/HB15) and 
KANADI TF (KAN1-KAN3) were shown to be involved in the 
secondary walled cell type formation and patterning in roots and 
stems (Baima et al, 2001; Emery et al, 2003; Kim et al, 2005; 
Ilegems et al., 2010). Three of the homeodomain TF mutants ana- 
lyzed in our study exhibited SCW phenotypes. The blh.6 mutant 
had less lignified SCW mainly in the xylary and interfascicular 
fibers, whereas in the hb5 mutant the fibers in both fascicular and 
interfascicular regions were heavily lignified. In the hbl5 mutant, 
both regions were also highly lignified but in addition ectopic 
lignification was observed in the parenchymatous cells adjacent 
to fiber and xylem cells (Figures 4B,D). This suggests that HB15 
represses the SCW formation program rather than only pro- 
mote the xylem cell differentiation as was concluded from earlier 
studies where down-regulation of CAN/HB15 stimulated xylem 
production, and over-expression of a miR166-resistant HB15 
(gain-of-function mutant) resulted in reduced xylem formation 
(Kim et al., 2005). Co-expression analyses revealed interesting 
clues for BLH6 which was co-expressed with genes involved in the 
biosynthesis of the three main polymers i.e., cellulose, xylan, and 
lignin as well as with the master regulator MYB46 and its clos- 
est paralog MYB83. Together with the hypolignified phenotype of 
the mutant and the thinner SCW particularly in the fibers, this 
further supports a role of BLH6 as an activator of the whole SCW 
program. 

AP2 ERF TF (AT3G14230) was identified in all four SCW- 
related transcriptomic datasets and exhibited high and preferen- 
tial expression in xylem, but the corresponding mutant had no 
detectable cell wall phenotype. Twelve members of the AP2 ERF 
TF family were highlighted by our in silico approach, seven of 
which had high and preferential expression in xylem and another 
(AT5G61590) was strongly induced during in vitro tracheary ele- 
ment formation. Although this family was the second most highly 
represented TF family just after the MYBs in the 80 candidate 
list, none of its members have yet been shown to be directly 
involved in the regulation of SCW formation. This family there- 
fore deserves more attention especially because it was recently 
reported that ethylene regulates cambium activity and promotes 
secondary walled xylem formation (Love et al., 2009). 



Some members of the auxin-dependent TFs Aux/IAA and 
ARF families have been shown to be involved in vascular tissue 
formation. For example, loss-of-function in ARF5/MP (Hardtke 
and Berleth, 1998) and gain-of-function in IAA12/BDL (Hamann 
et al., 2002) resulted in reduced and discontinuous vascular for- 
mation. These TF families were also highly represented within 
the 80 candidates with seven and two members for Aux/IAA and 
ARF, respectively. IAA9 was a very promising candidate found 
in the four transcriptomic datasets, highly and preferentially 
expressed in xylem and during tracheary elements differentia- 
tion. Unfortunately the corresponding mutant was unavailable 
at the time this work was performed. T-DNA insertion mutants 
corresponding to ARF4, ARF6 and 7AA28, and an IAA11 RNAi 
transgenic line were analyzed here but did not show any obvious 
SCW phenotype. This is very likely due to their functional redun- 
dancy as reported in previous studies (Okushima et al., 2005; 
Overvoorde et al., 2005). The creation of double/triple mutants 
of these paralog genes might be necessary to further assess their 
involvement in SCW formation. 

The hypolignified lines blh6 and zinc finger TF displayed earlier 
flowering time as compared to control whereas the hyperlignifi- 
cation line hb5 exhibited delayed flowering time. Two previous 
studies demonstrated that flowering induction time was deter- 
minant for xylem expansion and SCW formation in Arabidopsis 
hypocotyls and roots. Some major QTLs for SCW thickening dur- 
ing xylem expansion and fiber differentiation correlated tightly 
with a major flowering time QTL. In addition, transient induc- 
tion of flowering at the rosette stage promoted SCW thickening 
and xylem expansion (Sibout et al., 2008). Double mutant of 
two flowering time genes socl ful showed a synergistically delayed 
flowering time and a dramatically increased SCW formation with 
wood development present throughout all stems and to a much 
larger extent than any Arabidopsis mutant described to date 
(Melzer et al., 2008). Collectively these results suggest that the 
flowering induction is coupled with the SCW thickening program 
and xylem formation. 

In conclusion, we described here a post-genomic approach 
that enabled us to propose a list of 80 promising candidate genes 
potentially regulating SCW formation and/or lignification. Many 
of the available mutants analyzed did not provide any detectable 
SCW phenotype and complementary approaches (overexpres- 
sion, using different alleles, dominant repression, or multiple 
mutants) are now necessary to further characterize their func- 
tion. However, the six TFs of which mutants exhibited clear lignin 
phenotypes, further highlight the complexity of the regulatory 
network controlling SCW formation. Their in depth functional 
characterization should allow a better understanding of the regu- 
lation of lignification and SCW formation which may ultimately 
be used to improve the saccharification potential. 

MATERIALS AND METHODS 
CROSS-COMPARISON OF MICR0ARRAY DATASETS 

Four in house microarray datasets were generated in our labora- 
tory. In brief, datasets are from watl T-DNA Arabidopsis mutant 
CATMA microarray (Ranocha et al, 2010); EgMYBl (Legay 
et al., 2010), EgMYB2 over-expressed in Arabidopsis (unpub- 
lished), and orthologs of Eucalyptus xylem expressed genes 
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(Rengel et al., 2009). Publicly available microarray datasets were 
extracted from Genevestigator (https://www.genevestigator.com) 
(Hruz et al., 2008) by using Arabidopsis ATH1 22k array platform 
(7392 array datasets). 

PLANT MATERIAL AND GROWTH CONDITION 

The mutant lines were isolated from the T-DNA mutag- 
enized populations in the SALK collection (Alonso et al., 
2003) and from the RNAi transgenic plant populations in 
the Agrikola collection (http://www.agrikola.org). Seeds were 
obtained from the Nottingham Arabidopsis Stock Center (NASC) 
(http://arabidopsis.info/) and GABI (http://www.gabi-kat.de/). 
Homozygote lines were obtained from NASC or generated in lab 
and verified by PCR genotyping with gene specific primers and 
the respective left border primers of the T-DNA listed in supple- 
mentary Table SI 1. The transcript levels of each target gene in 
the six T-DNA insertion mutant were assessed (Figure SI) and 
the corresponding primers are listed in supplementary Table S12. 
Plants were grown in jiffy peat pellets then transferred to standard 
soil in culture room in short day conditions [9 h light, 200 u,mol 
photons m _1 s _1 , 22°C (day)/20°C (night), 70% RH]. The flow- 
ering time was considered from sowing day until the flower stem 
reached 20 cm in height. 

MICROSCOPY 

The histological comparative analysis of SCW between wild type 
and mutants was done at the stage of newly formed green 
siliques, about 2 weeks after bolting, when the inflorescence stems 
reach 20 cm in height. At this stage, the basal part of the inflo- 
rescence stem abundantly develops cells undergoing secondary 
wall thickening (xylem vessel cells, fascicular, and interfascicu- 
lar fiber cells). Lignin polymers are the characteristic components 
of SCW and are normally absent from primary cell wall, there- 
fore we used lignin deposition detection techniques to screen 
for SCW phenotype. Two methods were then chosen to detect 
the lignin polymers in the sections for microscopic observa- 
tion. Firstly we used the natural auto fluorescence of the aro- 
matic ring moieties on the subunits of the lignin polymer under 
UV-light exposition. Secondly, we used the phloroglucinol-HCl 
coloration which stains specifically lignin polymer precursors 
coniferaldehyde and p-coumaraldehyde in the SCW giving a red- 
purple color when observed under normal light. Cross sections 
of inflorescence stems at the basal end (100-150 u,m) were either 



observed using auto-fluorescence or stained with phloroglucinol- 
HCl. Auto-fluorescence was observed with a Leica microscope 
(excitation filter Bp 340-380 nm; suppression filter Lp 430 nm; 
http://leica.com). Phloroglucinol-HCl was directly applied on 
the slide. Images were recorded with a CCD camera (Photonic 
Science, http://www.pho tonic- science . co . uk) . 

CO-EXPRESSION ANALYSIS 

Three co-expression analysis tools were explored using 
Genevestigator (https://www.genevestigator.com), Arabidopsis 
co-expression data mining tools (http://www.arabidopsis.leeds. 
ac.uk/act/), and GeneCAT (http://genecat.mpg.de/). The results 
were presented using Genevestigator output tables and genes 
classified according to gene ontology semantic (Berardini et al., 
2004). We used Genevestigator Arabidopsis ATH 1 22k array plat- 
form with in absentia parameters that comprise all 7392 qualified 
datasets and is regardless of the underlying microarray datasets 
and the bait genes (i.e., all samples, condition-independent, and 
no-tissues specific bait genes), 50 was as "cut-off" threshold for 
co-expressed genes list. 
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